Dynamical magneto-electric coupling in helical magnets 
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Collective mode dynamics of the helical magnets coupled to electric polarization via spin-orbit 
interaction is studied theoretically. The soft modes associated with the ferroelectricity are not the 
transverse optical phonons, as expected from the Lyddane- Sachs- Teller relation, but are the spin 
waves hybridized with the electric polarization. This leads to the Drude-like dielectric function 
e(u>) in the limit of zero magnetic anisotropy. There are two more low-lying modes; phason of the 
spiral and rotation of helical plane along the polarization axis. The roles of these soft modes in the 
neutron scattering and antiferromagnetic resonance are revealed, and a novel experiment to detect 
the dynamical magneto-electric coupling is proposed. 
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The gigantic coupling between the magnetism and fer- 
roelectric properties is now an issue of keen interest 
Q&HSHEIzL&li A representative system of 
interest is i?MnC>3 with i?=Gd,Tb,Dy. For example, 
TbMnC>3 shows a ferroelectric moment P//c below a 
temperature Tfe = 28K(< T^eel — 42), and further- 
more changes its direction of electric polarization to a 
axis under the magnetic field H//b [Ell, 0, 0. Similar 
strong coupling behavior has been observed in N13V2O8 
|, i?Mn 2 5 (i?=Tb,Ho,Dy)[l[3, Bao. 5 Sri. 5 Zn 2 Fei 2 22 
8j . A common and essential feature of these compounds 
is that there are frustrations in the magnetic interac- 
tions. For i?MnC>3, Kimura et al. Q revealed that 
the increased GdFeC>3-type distortion of perovskite lat- 
tice leads to the further-neighbor exchange interactions 
and nontrivial magnetic structures. Later the neutron 
scattering experiment determined the spin structure of 
TbMnC>3; it shows the incommensurate collinear spin or- 
dering pointing along 6-direction for Tpe < T < T^eei, 
while the helical spin structure winding within the b — c 
plane occurs for T < Tfe with the helical wave vec- 
tor q/ /b [||. A similar helical structure is also observed 
in Ni3V 2 08 Q. These experiments point to key role of 
the non-collinear spin configurations such as helical spin 
structure, which are induced by frustrated exchange in- 
teractions, in producing the electric polarization and en- 
hanced magneto-electric coupling. 

A microscopic mechanism of the ferroelectricity of 
magnetic origin has recently been proposed by Katsura 
et al. which is based on the idea that spin cur- 

rent is induced between the noncolliear spins and hence 
is the electric moment due to Aharonov-Casher effect 
[IlHl2]| . This result can be regarded as an inverse effect of 
Dzyaloshinskii-Moriya interaction 1 13 . Phenomenologi- 
cal treatment of this mechanism |l 1[ and its extension 
to include electron-lattice interaction have been also 




FIG. 1: Schematic ground state configuration of spins(black 
arrows) and lattice displacements(blue arrows). 



reported. By now the static and/or ground state proper- 
ties of the magneto-electric coupled systems are well un- 
derstood. The obvious next step is to study their novel 
dynamic properties searching for the unique electromag- 
netic effects. In this context, we recall that the collective 
modes are a central issue in the research on ferroelectrics, 
where one of the transverse optical (TO) phonons softens 
toward the phase transition and is condensed according 
to the Lyddane-Sachs- Teller relation. This conventional 
view is not relevant in the case of multiferroics since the 
lattice displacement is not essential to the electronic po- 
larization. Therefore it is an important open problem to 
identify relevant collective mode that are responsible for 
ferroelectricity of magnetic origin. 

In this paper, we study exactly this question, i.e., the 
low energy dynamics of the system of the coupled spins 
and electric polarization. We found the new collective 
modes of spin and polarization waves that couple the di- 
electric and magnetic properties in a novel way. Experi- 
ments aimed at the detection of these collective modes in 
terms of the magnetic resonance, neutron scattering and 
the ac-dielectric measurement are also discussed. 

In ref. Katsura et al have shown that spin (su- 

per)current in noncollinear magnets j s (oc Si x Sj) leads 
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to the electric polarization P oc x j s with being 
the unit vector connecting two sites i and j . An effective 
Hamiltonian describing this coupling between the spins 
and the atomic displacement Uj, which represents the 
electric polarization including the displacement of elec- 
tronic cloud, is given as; 
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Here, the spin-lattice interaction A stems from the rel- 
ativistic spin-orbit interaction and corresponds to the 
DM interaction once the static displacement (Hj) is non- 
zero and breaks the inversion symmetry. Note that 
purely electronic polarization is possible even without the 
atomic displacement, but the latter is always accompa- 
nied with the inversion symmetry breaking. Therefore 
u m should be regarded as the lowest frequency repre- 
sentative coordinates relevant to the polarization, i.e., 
the transverse optical phonon, and the polarization p is 
given by e*u m with a Born charge e*. The term D(S V ) 2 
with positive D represents the easy-plane spin anisotropy. 
Here, we have assumed that the ground state spin con- 
figuration on the plane perpendicular to the helical wave 
vector is ferromagnetic and the low-energy excitations 
on this plane have a usual parabolic dispersion relation, 
i.e., oc \q±\ 2 . Hence we shall focus on the modes along 
the helical wave vector which is set parallel to the pr- 
axis shown in Fig. 1 . In H\ , the summation is taken over 
all the combinations of m, n and J{R m — Rn) denotes 
the Heisenberg interaction between S m and S n at R m 
and R n , respectively. First, we determine the classi- 
cal ground state configuration of spins and lattice dis- 
tortions. Since we consider the static configuration, the 
kinetic term in H3 can be negligible. The variation of H 
with respect to u m gives u m — ^e z x (S m x S m +i)- Hence- 
forth we assume that spins are on the easy plane, i.e., 
zx-plane and explicitly given by S£ = Scos(QRn + 0), 
SI = Ssia(QR n + 4>), SI = 0. By substituting this 
spin configuration into the energy e per spin is rewrit- 
ten as e = -S 2 J(Q) with J(q) = J(q) + ^-s'm 2 (qa). 
Here J(q) = J2 m J(R m )e~ lqRm is the Fourier trans- 
form of J(R m ) and a denotes the lattice constant. The 
wave number Q is determined to maximize J(q). Us- 
ing Q, the uniform lattice displacement is given as u m — 
— ^-sin(Qa)e x , and hence the static electric polarization 

is given as p = —pe x with p — e*u* n — e*^-sin(Qa). 
Now we consider the equations of motion for spins and 



displacements, and study the collective modes of our sys- 
tem. We introduce a rotating local coordinate system 
£, 77, £ such that the £-axis coincides with the equilibrium 
spin direction at each site, the £-axis is perpendicular 
to this direction in zz-plane and the 77-axis parallel to 
the y-axis (see Fig.l) |l6j. Assuming the small quan- 
tum/thermal fluctuations around the classical configura- 
tion, S*£ and can be regarded as a constant and the 
equation of motion can be written as S^ = S^H^ — SH%, 
S% = SHi - S£i?< where S = 5£ =const., and H n = 
H%, is the effective field acting on the n-th spin 
explicitly given by 



= 2S ( J{Q) + — sin 2 (Q a ; 



Hi = J2 2J(Rm - R n )Sicos(Q{R m - R n )) 

m 

+ ^sin 2 (Qa)(Si +1 +Si_ 1 ) 
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Here, we regard S^,S n and u v as sufficiently small and 
neglect their second- or higher order terms. We can also 
write down the equation of motion for u v . Introducing 
the Fourier components of S^, S^, and p y ml the equa- 
tions of motion are given by 



= -A{q)S\, 
S\ = B(q)S% - iXS 2 
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B(q) = 2S 



J{Q) - J(q) 
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We note that A(q) and B(q) satisfy the relations A(— q) = 
A(q),B(—q) = B(q), respectively, and A(Q) is equal to 
zero. 

From Eq.l JTj l. one can see the coupling between the spin 
wave modes and electric polarization. First S v and are 
canonical conjugate variables, and form a harmonic oscil- 
lator at each q in the rotated frame. This spin wave at q 
is coupled with the phonon u at q±Q, or u q is coupled to 
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S v /S^ at q±Q. Especially the uniform lattice displace- 
ment Uq is coupled to e~ l ^Sq — e^Stq (a =r],£), which 
corresponds to the rotation of both the spin plane and 
the direction of the polarization along the z-axis. This 
mode is the Goldstone boson with frequency lu = when 
D = 0, i.e., the symmetric case around z-axis. On the 
other hand, e^ 1 ^ Sq + e 1 ^ S"q corresponds to the rotation 
of spin plane along x-axis, which is decoupled to the po- 
larization but is gapped by the effective spin anisotropy 
introduced by the spin-lattice interaction. The spin wave 
mode at q = corresponds to the sliding mode, i.e., pha- 
son, of the spiral. This is decoupled from u and has zero 
energy at q = 0. 

Using Eq.JJJ) with canonical commutation relations 
[u q ,p q '] = iSq t - q > and [S^,S q ,] = iS5q.-qi, the matrix 
form of the retarded Green's function 



G R (AB:t-t') = -i6(t-t'){[A(t),B(t')]), 



G r (AB;lu) = 



2tt 
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(A,B = u q ,p q , Sq 1 ) is obtained. The imaginary part 
of the ac susceptibility is also defined as x"(AB;u>) = 
—ImG R (AB;oj). To be more explicit, we will discuss 
below some physically important cases. 

ac dielectric properties - The dynamical dielectric func- 
tion is given by e yy (iu) — 1 + 4:ir(e*) 2 G R (u uo; lu), where 
the e* is the Born charge corresponding to the displace- 
ment u v along y-axis, and 

,.,2_,.,2 

G r {u u q ;lu) 



2ttM(u; 4 - (wg + lu 2 )lu 2 + A{Q)Du%) ' 

(10) 

where lu p = A(Q)B(Q) is the frequency of the spin- 
plane rotation mode along x-axis, i.e., e~ l ^Sg + e^S^Q, 

and luq = \J k/M is that for the original phonon. This re- 
sponse function has the poles at ui± , which are explicitly 
given by 
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\ (ujI + cu 2 p ± ^4 + cu 2 ) 2 - 4A(Q)Dlu 2 ) (11) 



Assuming D,X <C lu 2 ,, one can estimate as oj_ = 
y/A{Q)D — VSSJD, and lu + = ujq. With this frequen- 
cies, one can write £ yy (cu) = 1 + J2± Ll >±I±/( uj2 ~ UJ '±) 
with the "oscillator strengths" I± being given by J_ = 



, and I + = 



2(e*)>^-^) 



Note that 



Mu-{uj\ - lu 2 _) ' + Muj + {u\-lu 
the integral — J °° dLue yy (uj) is given by (n/2)(I- + I + ). 
Here, the two modes contributing to the dielectric func- 
tion are (i) the phonon mode with the frequency tt>+ = too, 
and (ii) the z-axis rotation mode at w_ = -J A(Q)D. 
Usually loq is a high frequency and does not show any 
softening in the present case. The ferroelectricity is due 
to the spin ordering, which is hybridized with the polar- 
ization mode. In other words, the low frequency dielec- 
tric function is mostly due to the spin wave mode at ±Q. 



Note the low- frequency behavior of e yy (iS) is similar to 
the Drude form when D = 0, and the oscillator strength 
7_ (at lu = ui_) is enhanced as L_ ~ p 2 ^ J/D as D — > 0. 
This means that even though the spin-polarization cou- 
pling A and hence the static polarization p is small, the 
spin wave mode can contribute significantly when D is 
small. 

neutron scattering spectra - Next, we shall examine how 
the spin-polarization coupling affects the (inelastic) neu- 
tron scattering spectra. The intensity of the neutron 
scattered with the momentum transfer q is estimated as 
AmJ2 ae± G R {Sq*S« q ;uj), where the super- 



script a indicates the components perpendicular to q. In 
what follows, for the sake of simplicity, we set D = 0, <j> = 
0. In order to study the neutron scattering spectra, it is 
needed to rotate the spin coordinates back to the original 
laboratory system S a (a = x, y, z). We can easily derive 
the transformation formula given by 

G R (SqS : !L q ; lu) — - ^2 G R (S q+p S^ q _ p ,;uj), 

p=±Q 

p' = ±Q 

G R {Sy q Sl q] Lu)^G R {S^ q Sl q -uo), 
G R {S* q S*_ q ;u) = \ sgn(pp')G R (Sl +p Si q _ p ,;uj). 



p=±Q 
p' = ±Q 



(12) 

For example, at the helical wavevector q = (0, 0, Q), one 
can read the neutron scattering spectra by using the 
above transformation formula l|12|) and the list of pole 
positions and corresponding intensities summarized in 
Table |Ja) and (b). Note that the Green's function is 
written as G R (cu) = J2i a i/ ( cj2 — ^f) using the intensity 
dj and the pole position cjj. Now 4 modes are expected 
to contribute to the neutron spectrum |l7j : (i) phason 
mode at u) = 0, (ii) x-axis rotation mode at u) = lu p , (iii) 

phonon mode with q = at lu — \J lu 2 + uj 2 , and (iv) 
phonon mode with q = Q at loq. The former two, i.e., 
(i) and (ii), are of magnetic origin, and can be detected 
without the spin-polarization coupling, while the weights 
of the latter two, i.e., (iii) and (iv), are borrowed from the 
spin wave, and hence their relative magnitudes to the for- 
mer ones are roughly estimated as; (iii)/(i)~ {oj p /loq) 2 , 
and (iv)/(i) - X 2 /(k,J). 

antiferromagnetic resonance - We shall briefly discuss 
the antiferromagnetic resonance in our system. Here, 
D = 0, and = are also assumed. Magnetic resonance 
experiments pick up uniform magnetic excitations and 
corresponding Green's functions are G r (SqSq- u>), (a = 
x,y,z), where axis a corresponds to the direction of 
applying oscillating magnetic field. Using again the 
transformation formula Ijl2|l . one can see the pole 
of G r (SqSq ; lu) occurs at lu = lo p with the inten- 
sity SB(Q)/{Ait), while that of G r (S$S§;lo) at lu = 
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TABLE I: The pole positions u>i and corresponding intensities 
ai for each Green's function. 



y/^o + ^Jp with the intensity SB(Q)/(4n). Note here 

that there is no dynamics of Sq since it is a conserved 
quantity. The phonon mode with q = ±Q has a high 

frequency u) = y/u>Q + u)^ and should be difficult to be 

observed experimentally. The x-axis rotation mode of 
low frequency (u> = oj p ) is not coupled directly to the po- 
larization, while the resonance intensity is due to the spin 
anisotropy induced by the spin-lattice coupling A and/or 
the spin anisotropy D as discussed earlier [lsj . 

Now we turn to the realistic estimation of the param- 
eters in our theory extracted from the experiments on 
i?Mn0 3 . The exchange coupling constant can be es- 
timated from the spin wave dispersion observed in the 
neutron scattering experiment |19j . They estimated the 
exchange coupling J\ between the nearest neighbor spins 
as 85Ji 9 meV for PrMn0 3 and 85Ji 2.4 meV 
for TbMn0 3 . On the other hand, the anisotropy en- 
ergy, which corresponds to D in our theory, is around 
0.4meV for both compounds. From the ESR spectra 
of Jahn-Teller distorted LaMnO 3+5 (0 < 8 < 0.07), 
Dzyaloshinskii-Moriya coupling is estimated as T> ~ 
O.lmeV [i^ and hence we can evaluate the spin-lattice 
coupling A ~ lmeV/A ^^|. Assuming k ~ leV/A 2 
we can obtain the static displacement u m — 10~ 3 A. An- 
other important quantity is the Born charge e*. This 
can be estimated from the value of electric polarization 
with the above u m and the lattice constant a ~ 5A. 
The electric polarization along c-axis is determined as 
P c ~ 0.2/j,C/cm 2 for DyMn03 |3] and hence we can es- 
timate the Born charge e* as 16e where e is the bare 
unit charge. This result, i.e., the Born charge is much 
greater than the unit charge, indicates that the electric 
polarization mainly consists of the displacement of the 
electronic cloud. At this moment, let us check the consis- 
tency of our estimation with another recent experiment 
on dielectric response |2l|. They observed the peak of 
Ime at around 20cm~ 1 with the magnitude of 1 ~ 2 in 
GdMn0 3 and TbMn0 3 . This 20cm- 1 is identified with 
o;_ in Eq.(12), and the integration of — \m£ vy (ui) over w 
gives the "oscillator strength" (jr/2)I-. From their re- 
sults, we can estimate i_ ~ 12cm _1 . On the other hand, 
w_ and J_ are independently estimated as cj_ ~ 10cm _1 



and I- ^4cm _1 from the parameters determined above. 
Here, we have assumed that S — 2 and Qa = tt/4. Both 
uj- 's and J_ 's, are of the same order and hence we can 
conclude that our phenomenological model well describes 
the magneto-electric coupling in i?Mn0 3 . It turned out 
that the coupling between the spin wave and the trans- 
verse phonon is weak, and is characterized by the ratio 
A 2 /(kJ) ~ 10~ 3 . Therefore, it would be rather difficult 
to see this coupling in the neutron scattering experiment, 
i.e., the magnetic scattering by the phonon modes. How- 
ever, note that the low energy dielectric response reflects 
the spin wave mode as the Goldstone boson, i.e., z-axis 
spin rotation mode, and it would be enhanced more if 
the spin anisotropy is weaker. 

Lastly, we propose a novel experiment to see directly 
the dynamical magneto-electric coupling. Namely we can 
excite the spin wave by electric field. These responses 
are described by the off-diagonal element of the Green's 
function; 

G r {uIS*-u) = -G r {S*uI-uj) 

1 i\S 2 s\n{Qa)uj 

~ 2ttM uj* - (wg + w2) w 2 + A(Q)DujI ' ( ' 

Proposed experiment requires to put the sample between 
the two electrodes and apply the a.c. electric field E y of 
frequency uj along y-direction. Then the magnetization of 
the same frequency u can be detected at the edges of the 
sample polarized along z-direction. Now oj_ = 20cm _1 
corresponds to the Tera-Hertz region, and the highest 
intensity of the electric field there is ~ 10 3 V/cm. Putting 
this value and all the other parameters estimated above 
with overdamping 7 ~ c<j_ = 20cm, , we obtain the 
estimation for the magnitude of the a.c. magnetization 
as m z ~ 0.4 x 10 _4 /ZB/Mn atom, which can be detected 
by Kerr rotation spectroscopy. We should note here that 
the magnitude of m 2 is much enhanced when 7 is small. 
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